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ABSTRACT 

— M e thod A method of generating a hybrid grid allowing modelling of a 
heterogeneous formation crossed by one or more pipes such as, for example, an 
underground formation where one or more wells have been drilled, in order to 
form a representative models for example of fluid flows in th i s m e d i um the 
formation in accordance with a defined numerical pattern is disclosed . 

The method e ss e nt i a ll y comprises associating a first structured grid for 

gridding of the heterogeneous medium respect i ng tho regarding discontinuities 
thereof with a second structured, radial type grid for gridding of a zone around 
each pipe or well, which allows to-better rospoot particu l ar constraints linked with 
flows in tbis4he_zone, and transition of non-structured grids that are interposed 
between the first grid and each second well grid. Various grids are combined, 
each with its own formation, representation and exploration methods, structured 
jgrids which are advantageous i n that thev fac il itate in facilitating control and 
comprehension of the reservoir images formed and more flexible non-structured 
grids for gridding of complex zones. An application is simulation of 
— Appl i cants: hydrocarbon r e s e rvoir simu l ators for e xamp le reservoirs . 
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BACKGROUND OF THE INVENTION 
F I ELD OF THE INVENT I ON Field of the Invention 
The present invention relates to a method of generating a hybrid grid 
allowing modelling of a heterogeneous formation crossed by one or more pipes. 
5 The method of more particularly applied to formation of a grid suited to an 

underground reservoir crossed by one or more wells, in order to model 
displacements of fluids such as hydrocarbons. 



BACKGROUND OF THE I NVENTION Description of the Prior Art 
1 0 Grid generation is a crucial element for the new generation of reservoir 

simulators. Grids allow to d e scr i b e description of the geometry of the geologic 
structure studied by means of a representation in discrete elements wherein 
simulation is performed according to a suitable numerical pattern. Better 
comprehension of physical phenomena requires 3D simulation of the multiphase 
1 5 flows in increasingly complex geologic structures, in the vicinity of s e v e ral types 
of singularities such as stratifications, faults, pinchouts, channels and complex 
wells. All this complexity has to be taken into account first by the grid which has 
to reproduce as accurately as possible to geologic information in its 
heterogeneous nature. 
20 Grid modelling has made great advances during the past few years in 

other fields such as aeronautics, combustion in engines, structure mechanics, 
etc. However, the gridding techniques used in the other fields cannot be applied A 
as they are A to the petroleum sphere because the professional constraints are not 
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the same. For example, in reservoir simulation, the numerical patterns are 
constructed from control volumes in order to better respect the mass 
conservation in the case of transport equations of hyperbolic nature. The grid 
must be a « b l ock c e nt e r e d» "block-centered" type grid, irer that is the nodes 
5 must be inside each layer and the boundaries of each block must follow the 

interface between the layers. Now, if this constraint was not taken into account, 
the nodes would naturally be placed along the faults and along the stratification 
boundaries. The consequence of this would be that these interfaces would pass 
through the control volume that is used. The saturation, constant in the control 

10 volume, could not r o spoot th e n th e consider discontinuity and the results would 
not be accurate. It is therefore necessary to develop new techniques that are 
better suited to tho potro l oum sph e r e requirements of the petroleum field . 

Cartesian grids, which are commonly used in current commercial 
simulators, are unsu i t e d not suitable for solving these new problems posed by 

15 the development of petroleum reservoirs. Cartesian grids, based on 
parallelepipedic elements, do not allow representation of such complex 
geometries. 

There is a well-known method of generating structured 3D hexahedral 
grids of CPG (Corner-Point-Geometry) type which r e sp e cts considers t he 
20 geometry of the bodies, It is described in pat e nt FR - French patent 2,747,490 
(US -patent 5,884,564) filed by the app l ic a nt assignee and also in the following 
publication: 



Bennis Ch. Et al. «On e Mor e St e p i n Gocad Stratigraph i o Gr i d 
G e n e ration^ "One More Step in Gocad Stratigraphic Grid Generation : Taking 
into Account Faults and Pinchouts" ; SPE 35526, Stavanger, 1996. 

This grid type is more flexible than Cartesian grids because it consists of 
5 any hexahedral elements that can be degenerated. It strictly respects the 
horizons, the faults and it allows to r e pr e s e nt representation of certain 
unconformit ie s nonconformities such as pinchouts because its construction is 
based on these elements. However, this type of grid does not allow to so l v e 
solution of all the geometric complexities such as, for example, circular radial 
10 grids around complex wells. It is possible to form separately the grid of the 
reservoir and the girds around the wells but it is difficult to represent several 
objects in the same CPG type reservoir grid because of connection problems 
linked with the structured nature of the grid. 

Another approach is also known where 3D grids only based on tetrahedral 
15 Delaunay elements, with a circular radial refinement around the wells, are 
automatically generated. The advantages of such an approach is that it is 
entirely automatic and does not practically not require the user's attention. 
However, this method has drawbacks which make the results obtained difficult to 
use: 

20 - there are on average five times as many grid cells as in a CPG type grid 
for the same structure, which is very disadvantageous for simulation calculations, 
unlike the structured grids which are easy to visualize, to explore from the 
inside and to locally modify interactively, it is very difficult and sometimes 



impossible to properly control the tetrahedral grids because of their size and 
especially because of their non-structured nature. This poses problems for 
validating the grid from a geometric point of view as well as for understanding 
and validating the result of a simulation of this type of grid. 
5 Other approaches are also well-known, which allow to g e n e rat e 

generation of grids, notably grids based on control volumes generated from a 
triangulation (VoronoT and CVFE), associated with techniques of aggregation of 
the triangles (or tetrahedrons) into quadrangles allowing the number of grid cells 
to be reduced. Although promising results were obtained with these new grids, 

10 precise representation of the geologic complexity of reservoirs and wells remains 
a subject for research and development. In fact, these grids are rather 2.5D (i.e. 
vertically projected) and their 3D extension appears to be very difficult. Despite 
their hybrid aspect, they remain entirely unstructured and would therefore be very 
difficult to manage and to handle in real 3D. Furthermore, taking account of real 

15 3D faults and deviated wells would greatly increase this difficulty. 



SUMMARY OF THE INVENTION 
The method according to the invention allows to g e n e rat e generation of a 
hybrid grid suit e d to applicable to a heterogeneous formation crossed by at least 
20 one pipe of known geometry (such as an underground reservoir crossed by one 
or more wells), in order to form a model representative of fluids flows in tWs-the 
formation in accordance with a defined numerical pattern, the structure of the 
formation being known a priori from available data acquired by means of in-situ 



measurements, analyses and/or interpretations of formation images (seismic 
images for example, in the case of a reservoir). 

The method i s charact e r i z e d i n that i t comprises associating a first 
structured grid for gridding the formation wh ile r e sp e cting th e with consideration 
5 of discontinuities thereof with second structured, radial type grids for gridding of 
the zones around the wells, these second grids allowing to r e sp e ct consideration 
of constraints linked with the-flows in the wells, and non-structured transition 
grids between the first grid associated with the formation and the second grids 
associated with the wells. 

1 0 Gridding of the heterogeneous medium is obtained A for example^ by 

importing each second structured grid into a cavity formed in the first structured 
grid, the size of this -the cavity being sufficient to allow formation of a non- 
structured transition grid between the first structured grid associated with the 
formation and the second structured grid associated with each well. 

1 5 The non-structured transition grid based on any polyhedrons and 

canonical polyhedrons such as pentahedrons, tetrahedrons, pyramids, etc., can 
be formed by r e sp e ct i ng consideration of constraints linked with the numerical 
pattern. 

The non-structured transition grids are advantageously modelled with the 
20 structured well grids by applying a forma li sm technigue known in the art, referred 
to as the _*:<g e n e ral i z e d map» forma li sm "generalized map" technigue . the grid 
of the formation being structured metrically, globally or in faulted blocks. 
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The global hybrid grid is thus obtained by combination of several grid 
types : a structured reservoir grid, a radial grid around each well, also structured, 
and non-structured transition grids which connect the previous two grid types. 
Each one of these grids has its own formation, representation and exploration 
5 methods. The structured aspect is thus degraded only at the points where this is 
strictly necessary. This «obi e ct» "object" approach affords both the advantage 
of structured grids for control and comprehension of the reservoir and the 
flexibility of non-structured grids in complex zones. Complexity is introduced only 
where it is strictly necessary. The independence of these gridding modes 
10 therefore allows separate extraction, management and representation of the well 
grids and the interstitial grids includes in the reservoir grids. 

Using a reservoir simulator of a well-known type, such as ATHOS™ or 
SCORE™ for example, for a reservoir provided with a hybrid grid obtained by 
means of the method, allows production simulations to be performed. 

15 

BRIEF DESCRIPTION OF THE DRAWINGS 
Other features and advantages of the method according to the invention 
will be clear from reading the description hereafter of non limitative examples, 
with reference to the accompanying drawings wherein: 
20 - F i gur e Fig. 1 shows a diagrammatic example of a hybrid grid of a reservoir 
crossed by two wells, consisting of a first structured grid for the reservoir, a 
second structured grid for the zones around the wells and transition grids 
between the two grid types T i 



• f 

F i gur e Fig. 2 shows an example of a structured grid of a faulted reservoir T ; 

F i gur e Fig. 3 shows an example of a radial grid around a vertical well T ; 

F i gur e Fig. 4 shows an example of grid of a horizontal well T ; 

F4gwe -Fig. 5 shows, in 2.5D, an example of a gridded reservoir where 
5 cavities are provided for gridded wells, before the stage of creation of non- 
structured interstitial grids intended to connect them together?; 

Rowe -Fig. 6 shows five wells provided each with a radial grid, included in 
a gridded reservoir, by means of non-structured transition grids based on any 
polyhedral grid cells T ; 
10 - F i gur e Fig. 7 shows separately the structured reservoir grid with the 
cavities provided therein in order to include the additional elements : gridded 
wells and interstitial grids T ; 

F i gur e s Figs. 8A, 8B, 8C, 8D show elementary well grids v i sua li z e d 
illustrated individually according to different modes T ; 
15 - F i gur e s Figs. 9A, 9B, 9C, 9D show various elementary transition grids, 

visual i z e d illustrated individually according to different modes, also, allowing their 
integration into the reservoir grid T ; 

ftgwe -Fig. 10 shows an example of a hybrid reservoir grid with transition 
grids consisting for example of pentahedrons, between the reservoir grid and 
20 several well grids T ; 

F i gur e Fig. 1 1 shows a model for a matrical representation of a structured 
grid around a well T i 
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Rawe- Fig. 12 is a graphic representation of a connection between 
strands, within the scope of the modeling technique referred to as the 
«q e n e ra li z e d maps» "generalized map" technique used for generating non- 
structured grids 7 ; 

5 - Figur e s Figs. 13A t 13B are graphic representations of connections by 
means of simple arcs T ; and 

Figur e s Figs. 14A, 14B are graphic representations of connections by 
means of double or triple arcs respectively^ 



10 DETA I LED DESCR I PT I ON DETAILED DESCRIPTION OF THE INVENTION 
Modelling of the reservoir is obtained by combining elementary grids of 
different types. Each grid is considered to be a full object with its own data 
model, its own generation methods and its own representation methods. 
Generation is carried out in stages with addition/subtraction of grids. 

15 1 ) In order to represent the reservoir as a whole, an i, j, k structured grid 

of a known type known to sp e c i al i sts , referred to as CPG (Corner Point 
Geometry), as described in the aforementioned pat e nt FR - French patent 
2,747,490 for example is used for e xampl e. The reservoir can be faulted with 
downcreep of a block in relation to the other. The major horizons and faults are 

20 first modelled by continuous surfaces from data resulting from an interpretation of 
seismic images of the subsoil or from data obtained during drilling (well markers). 
The geologic structure is then divided into faulted blocks resting on these 
surfaces. These blocks are individually gridded, then reassembled. Gridding of 
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a block first consists in gridding the edge surfaces, then the inside is populated 
by transfinite interpolation of the edge surface grids. Relaxation techniques are 
then applied to the edge surfaces and to the inside so as to harmonize and the 
regulate the grid. The grid thus obtained strictly r e sp e cts considers the 
5 horizons, the faults and allows to r e pr e s e nt c e rtain unconformiti e s representation 
of nonconformities such as pinchouts. It meets all the constraints of geologic 
nature. 

2) A well trajectory is drawn synthetically or imported. A structured radial 
grid is then generated around each well in order to take into account of th e 

10 particular constraints linked with the flows in the vicinity of these- the wells. 

In the example shown in Fig. 3, the structure grid around a vertical well is 
of circular radial type. It is also a CPG type grid. Its generation first consists in 
sampling a disc at r, 0 in the horizontal plane. The 2D grid thus obtained is then 
projected vertically upon the various layers of the reservoir grid. Here, the i, j, k 

15 of the matrical structure correspond to the samplings at r, 0 and z respectively 
(see Fig. 11). 

The grid around a horizontal well (Fig. 4) is i, j, k structured, it is of the 
same type as that of the reservoir, except that a well cannot be faulted. It is also 
obtained by projecting vertically upon the various layers of the reservoir grid a 2D 
20 grid belonging to a horizontal plane. 

3) This radial grid is the inserted around the-oi^around each well in the 
global reservoir grid. A cavity is therefore first created in the reservoir grid by 
deactivating all the grid cells in contact with well grid cells (Figs. 5, 6). The space 
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freed between the reservoir grid and the well grid must be sufficient to allow 
convenient formation of a transition grid. It can represent for example about the 
equivalent of two grid cell layers. 

4) A non-structured transition grid is then generated in this cavity (Fig. 7) 
5 in order to connect the structured radial grid around the well to that of the 

reservoir best respecting the constraints linked with the numerical pattern. The 
user can deactivate the grid of a well any time by reactivating the grid cells of the 
corresponding cavity in the reservoir grid. 

The transition grid can for example consist of polyhedrons with any 
10 number of sides or canonical polyhedrons (tetrahedrons, pentehedrons, 
pyramids, etc.) according to the numerical pattern used, without the overall 
hybrid approach proposed being affected. 

Example of modelling of a hybrid grid 
15 The reservoir grid and each well grid are modelled, for each faulted block 

of the reservoir, by matrical structures of points or cells each comprising e ach 
eight points. Because of the structured nature of the grids, the topological links 
between the various grid cells are implicitly contained in the matrical structure. 

Transition grids are more difficult to manage because of their non- 
20 structured nature and because they can contain polyhedral grid cells whose 
number of sides varies from one cell to the other. 

An advantageous solution for facilitating management of this new grid 
type, allowing to browse it and to surf it efficiently, cons i sts i n us i ng usesthe 
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topological model referred to as the «g e n e ra li z e d maps or G maps» 
"generalized maps or G-maps" . This model known to specialists i s and for 
example described by: 

Edmont J.: «_7\ Combinatorial Representation for Polyhedral 
5 Surfaces^!, Notice Amer. Math. Soc, 7, 1960, or by: 

Fortune S., 1992: Voronoi d i agrams Diagrams and Delaunay 
tr i angulat i ons Triangulations . pp. 225-265 of D.Z. Du & F.K. Hwang (eds.), 
Computing in Euclidean Geometry, 2 nd edn. Lecture Notes Series on Computing, 
vol. 4, Singapore, World Scientific. 
10 Generalized maps are based on a formal algebraic approach that is briefly 

r e m i nd e d described hereafter. 

In 3D, the elements which constitute a generalized map are (D, a 0 , ai, 02, 
a 3 ), where D is a finite set of elements called strands, and elements {a} are 
involution on D type functions, associating the strands two by two at most, which 
1 5 are therefore conveniently referred to as links. Figs. 1 1 to 14 show concr e t e 
specific geometric representation examples. Link a G is in the form of a dotted 
segment (Fig. 11) and links ch, a 2 and a 3 in the form of arcs, respectively simple 
(Fig. 12), double (Fig. 13) and triple (Fig. 14). 

According to another known approach, generalized maps are considered 
20 as graphs whose strands from the nodes and the links form the arcs : link a 0 

between two strands can be used for representing the edge of a side, links ai for 
connecting two edges of a side, links a 2 for linking two sides of a cell together 
and links a 3 for stick i ng connecting two cells together. 
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This -The model of generalized maps involves a small number of formal 
objects and an operation which, by associating additional information with a 
topology, allows to l ocat e location of the objects defined in space and to account 
accounting f or their appearance, which is referred to as "plunge^ and, in the 
5 present case, "plunge in a 3D space". 

It affords the advantage of being independent of the dimension of the 
objects. All the objects can be represented with the same data structure and 
handled with the same methods. This approach makes it possible to handle 
objects created with heterogeneous topological models. It is therefore well- 
10 suited for implementing the method according to the invention with i ts stag e of 
creation of a non-structured grid linking two different structured grids together. 

The generalized map conc e pt technique f or modeling the transition grid is 
applied by creating a c e rtain number of objects of different types which refer to 
one another. These objects mat e ria li z e define the topological network and its 
15 various plunges in a 3D space. Coh e r e nt l v Specificallv . in order to allow browsing 
the grid, a topological network is constructed parallel to the geometric data 
commonly handled in a grid, the points, the edges, the sides and the cells. 
Furthermore, crossed links are established between the topological network and 
geometric data. 

20 

Objects 

The various objects handles within the scope of the application performed 
here of the generalized map method are as follows: 
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1 ) The Transition Grid object which contains all the topology, the 
geometry and the physical data. It cons i sts of is_a GMap type object which 
represents the topological network and of a Plunge type object which 
mat e ria li z e s defines the plunge in the physical world according to ou^-the 
5 application. 

2) The GMap object 

The topological model is entirely contained in a graph cons i sting of a list of 
Strands connected to one another. Any operation performed on the generalized 
10 map amounts to an operation of the Strands network. The GMap object type has 
methods allowing easy circulation in the topological network representing the 
grid, irer that is to go from one Strand to another. 



3) The Strand object 

15 Each Strand is defined by four references to other Strands (corresponding 

to links a 0 , ch and a 3 ) and by four other references to the plunge in the 3D space, 
notably at a Point, an Edge, a Side and at a Cell to which it is connected. 



4) The Plunge object 
20 It is defined by four lists : 

a list of Points (the grid points), it is the plunge of dimension 0 of the 

GMap, 

a list of Edges (plunge of dimension 1 of the GMap), 
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a list of Sides (plunge of dimension 2 of the GMap), and 
a list of Cells (plunge of dimension 3 of the GMap). 

The Plunge object also contains its own methods of creating and handling 
the data it contains according to the use that is made thereof. Boo i dos. tho The 
5 GMap is created from its Plunge. 

5) The Point object 

A Point is defined by its coordinates x, y, z and by a list of attributes, 
notably scalar or real petrophysical values that are associated therewith. 

6) The Edge object 

It is defined by a reference in the GMap to a Strand which represents an 
end of the Edge. This gives a preferential access to the topological network and 
simultaneously allows te-ga -going from the plunge to the strands graph. For 
example, link a G of the Strand in question leas to the Strand representing the 
other end of the Edge. It is furthermore defined by a list of attributes, notably 
scalar or real petrophysical values that are associated therewith. 

7) The Side object 

20 This type of object allows to handl e handling directly the interfaces 

between the cells as well as the outer sides of the grid. A side is defined by a 
reference in the GMap to a Strand which represents a vertex of the edge polygon 
of the Side. This also gives a preferential access to the topological network. The 
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Strands representing the other vertices of the polygon are accessible by iterative 
applications of the relation a 0 och by starting from the initial Strand and eventually 
coming back to this same Strand. It is also defined by a list of attributes relative 
to the Side (for example, scalar or vector physical values). 

5 

8) The Cell object 

The type of Cell object is defined by a list of references in the GMap to 
Strands, each on representing a half-Side of the Cell. This allows access to the 
topological network from the Cells. It is also defined by a site (coordinates of the 
10 center of mass of the Cell) and by a list of attributes specific to the Cell or to its 
site (scalar or vector petrophysical values for example). 



Graphical representation and exploration 

A graphical representation is a very efficient and even essential means for 
1 5 controlling and validating the construction of a grid and the simulation results. 
Concerning the constructions, the geometry of the grid generated is generally 
first visually controlled. If this is not sufficient, local or global quality criteria with 
which statistics are established can be calculated and visualized on the grid by 
means of a eeJom ^color scale. Flow simulation cons i sts i n ca l cu l ating calculates 
20 the variations with time of certain petrophysical parameters by taking aocount of 
accounting for the hypotheses that initially condition the flows. Simulation 
validation also involves visualization of these parameters on the grid (preferably 
by means of a colour scale). As the involved grids concornod are 3D grids, tools 
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allowing to e xp l or e exploration of the grid from the inside by visual browsing are 
required. Graphical representation and browsing in the grids, presented 
hereafter, are a good illustration of the flexibility and the modularity of the hybrid 
approach proposed and of the efficiency and the adequacy of the data model 
5 selected. 

The hybrid grid T Js considered as a set of independent entities^ the -The 
elementary grids T is constantly v i sua li z e d seen in a main window. The user can 
select at any time an elementary grid and v i sual i z e see it with its specific 
methods in a secondary window which contains only the selected elementary 

10 grid se l ect e d . Actions on the elementary grid have automatically repercussions 
on the entire hybrid grid v i sua li z e d seen in the main window. An elementary grid 
can thus be v i sualiz e d seen and explored as a full entity and it can be viewed in 
the global context. The visualization methods differ according to whether the 
elementary grid is structured (reservoir grid and well grids) or not (transition 

15 grids). 

Examples of functionalities specific to structured grids 

In the case of a structured grid (reservoir and well), visualization is simple 

and conventional. It consists in two main functionalities: 
20 - visualization of the external envelope of the grid with the possibility of 

peeling it in the 3 directions i, j, j separately, 

simultaneous or separate v i sua li zation viewing of three matrical cell slices 

i=cste, j=cste and k=cste, with the possibility of moving them -slices in the block. 
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Examples of functionalities specific to non-structured grids 

In the case of non-structured transition grids, other, more elaborate 
visualization modes are preferably used. Five functionalities are mainly used: 
5 - visua li zat i on viewing of the external envelope with the possibility of 
concentric peeling, topological^ speaking, 

v i sua li zat i on viewing of the cells crossed by cutting plane orthogonal to an 
axis of coordinates x, y, z or any axis, 

visualizat i on viewing of the trace of the cells on the cutting plane, 
10 - v i sua li zat i on viewing of the grid sites when they are initially given and 

visu ali zat i on viewing of the grid cells in full or scattered mode. 

Of course, for the two grid types, it is possible to v i sua l iz e view a property 
or scalar value by means of a co l our color scale. 

All these functionalities require easy and optimum viewing of a non- 
15 structured grid. This is possible by using the forma li sm r e f e rr e d to as 
gonoral i zod map forma li sm generalized map technigue . 

Figs. 7 to 10 clearly illustrate the potential afforded by the hybrid grid 
m e thod propos e d, i . e . that is harmonious integration of a structured grid following 
a topological model (well grid) into another structured grid (reservoir grid) 
20 following a different topological model, by means of a non-structured transition 
grid. The independence of these models therefore allows extraction and 
separate representation of the well grids and of the interstitial grids included in 
the reservoir grid in order to represent, handle and explore this type of data. 



